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Abstract 

The emission from a radiating source embedded in a photonic lattice is calculated. The analysis 
considers the photonic lattice and free space as a combined system. Furthermore, the radiating 
source and electromagnetic field are quantized. Results show the deviation of the photonic lattice 
spectrum from the blackbody distribution, with intracavity emission suppressed at certain frequen- 
cies and enhanced at others. In the presence of rapid population relaxation, where the photonic 
lattice and blackbody populations are described by the same equilibrium distribution, it is found 
that the enhancement does not result in output intensity exceeding that of the blackbody at the 
same frequency. However, for slow population relaxation, the photonic lattice population has a 
greater tendency to deviate from thermal equilibrium, resulting in output intensities exceeding 
those of the blackbody, even for identically pumped structures. 

PACS numbers: 
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I. INTRODUCTION 



One of the many novel optical phenomena exhibited by photonic lattices is the modi- 
fication of spontaneous emission properties. ^'^ A photonic lattice can funnel radiation into 
narrow energy bands, where exceedingly high intensities at photonic lattice bandedges have 
been predicted theoretically and observed experimentally.^"^ A question is whether the peak 
intensities exceed those of a blackbody under similar experimental conditions. ^"^^ The an- 
swer is important for scientific understanding and can impact the development of new light 
sources. 

It is generally agreed that the higher photonic-lattice density of states will increase in- 
tracavity intensity. The debate concerns the output intensity in comparison with that of a 
blackbody. Arriving at an answer is difficult experimentally because it is difficult to ensure 
that the comparison is made under similar conditions. Theoretically, the difficulty lies with 
the treatment of the matter and optical aspects of the problem. The derivation of the mat- 
ter equations requires knowledge of the normal modes of the optical structure, preferably in 
the form of an orthonormal basis. However, such a basis set is not rigorously defined for a 
finite photonic lattice with outcoupling loss. This problem occurs also in laser theory, where 
one usually begins with the Fox-Li modes for a Fabry-Perot cavity with perfectly refiecting 
mirrors, and introduces a loss mechanism to represent the outcoupling.^^ Such a phenomeno- 
logical approach is inadequate for the present problem because of the inconsistency arising 
from separating the treatments of the eigenmode problem and the outcoupling cfi'ccts. 

The approach taken in this paper considers the photonic lattice and free space outside of 
the photonic lattice as one combined system (see Fig. 1). We follow the method of an earlier 
paper on the linewidth of a Fabry-Periot laser^^ in representing free space by a very large 
cavity. The photonic lattice is approximated by a series of semitransparent interfaces. We 
begin with discussing the one-dimensional geometry, ^'-"^^ which we will show to contain the 
essential fcatm^cs necessary for addressing our question. Section II discusses the equations 
and boundary conditions obeyed by the modes of our 'universe'. The determination of the 
eigenfrequencies and eigenfunctions requires the simultaneous diagonalization of a usually 
large matrix and the solution of a transcendental equation. A numerical procedure for a 
photonic lattice of arbitrary size and interface transmission is presented in Appendix A. 

In Sec. Ill, the radiation field is expanded in term of these large number of modes 
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and quantized. The radiation source is also treated quantum mechanically, as an inho- 
mogeneously broadened ensemble of two-level atoms confined within the photonic lattice 
structure. The equations of motion for the photon number and atomic populations are de- 
rived in this section. We choose a fully quantized (i.e., quantized matter and field) treatment 
based on Einstein's derivation of the Planck radiation law, which showed the importance 
of a consistent treatment of stimulated and spontaneous emission processes. Einstein was 
able to circumvent a fully quantized theory by using Wien displacement law, which applied 
only to emission in free space. For the photonic lattice, such a general relation does not 
exist. 

There arc several recent calculations of photonic-lattice emission where the emitting 
source is a classically described current. ^^'^^ An advantage of our treatment over these 
classical ones is that by paralleling the Planck radiation law derivation, our comparison of 
photonic-lattice and blackbody emission spectra is appreciably more straightforward. More- 
over, the fiuctuation-dissipation theorem^^ which is an essential assumption in the classical 
calculations, -"^^'^^ appears as a result in a fully quantized treatment because spontaneous 
emission is treated from first principles. 

Section IV describes a spectrometer model used to determine the emission spectrum. 
Section V uses the theory developed in the earlier sections to investigate emission from an 
active photonic-lattice that is excited by an external pump and allowed to equilibrate with 
a thermal bath via collisions. The radiation field spectra measured inside and outside the 
photonic lattice are described. Comparison of photonic lattice and blackbody emissions is 
discussed for equilibrium and nonequilibrium situations. 

Section VI summarizes the extension to a 3-dimensional geometry. The comparison 
between photonic lattice and blackbody spectra is made assuming a spherically symmetric 
photonic-lattice dielectric function. The 3-d treatment is important for three reasons. First, 
it proves that the theory can retrieve Planck's blackbody distribution in the absence of a 
photonic lattice. Second, it verifies the 1-d treatment in terms of containing the necessary 
physics for answering the question of photonic lattice versus blackbody thermal emission. 
Lastly, it points out the substantial increase in numerical demands with increase dimensional- 
ity, thus justifying our concentration on the 1-d analysis to facilitate physical understanding 
and tractability of numerics. 
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II. MODES OF THE COMBINED PHOTONIC-LATTICE AND FREE-SPACE 



SYSTEM 

In this section, the eigenmodes for a photonic lattice coupled to the outside world are 
derived using the model depicted in Fig. 1. The universe, which embeds the photonic lattice, 
is represented by a very large cavity with perfectly reflecting walls at 2; = and z — L. (End 
results are extrapolated by taking the limit L — > 00.) The photonic lattice is modeled as a 
series of coupled resonators with semitransparent interfaces. Following Spencer and Lamb,^^ 
the semitransparent interfaces arc described as very thin surfaces with very large dielectric 
constants. As an idealization, we use dielectric 'bumps' giving a dielectric permittivity 



e {z) = 60 



pi 



(1) 



where rj = 2-^ (1 — Tpi) /Tpi, Tpi is an effective transmission at each interface located at zj, k 
is the average magnitude of the electromagnetic field wave vector and Npi is the number of 
periods making up the photonic lattice. For brevity, we assume the background permittivity 
inside the photonic lattice to be that of vacuum eo. 

Using the above dielectric function in Maxwell equations gives the following differential 
equation for the eigenmodes of the combined photonic-lattice and free-space system: 

—life {z) = -HqE {z) nluk (z) (2) 

where ^0 is the permeability in vacuum, Q^. is the eigenfrequency and k labels the cigcnmode. 
The boundary conditions are obtained by first noting that the system is bounded by totally 
reflecting surfaces, so that 

UkiO) = UkiL) = (3) 
Integrating Maxwell's equations across the bump gives the boundary conditions, 

Uk{z^) = Uk{zj) (4) 

(4) - {zj) = -7)kuk{zj) , (5) 

where the superscripts — and -|- indicate the positions immediately before and after an 
interface, respectively. Integrating by parts (2) gives the orthogonality relation 

/ dz e {z) Uk {z) ui {z) = eo(5fe,; (6) 
Jo 



Plotted in Fig. 2 are examples of eigenfunctions for a six-period {Npi — 6) photonic 
lattice with effective interface transmission Tpi — 0.1. Most of the solutions are not resonant 
with the photonic lattice, so that mode amplitude is negligible inside the photonic lattice, 
as shown in Fig. 2 (a). Figures 2 (b) and 2 (c) show examples of photonic lattice modes, 
where the latter figure clearly depicts the first derivative discontinuities at the interfaces. 

We show in Figs. 3 and 4 that the model can reproduce the photonic-lattice properties 
relevant to our study. Figure 3 illustrates the formation of bands and bandgaps, by plotting 
the frequencies of the photonic lattice modes [i.e., modes depicted Figs. 2 (b) and 2 (c)] 
versus the interface transmission. Not plotted are the large number of free-space modes [Fig. 
2 (a)], with mode separation Q = rmrc/L — > as the system length L — > oo. At Tpi = 0, 
the photonic-lattice modes arc simply the modes of six uncoupled resonators, each of length 
a, i.e., they are A^p^-fold degenerate and have frequencies Q = mire/ a where m is an integer. 
The degeneracy is removed with couphng among sections of the photonic lattice. The result 
is groupings of states separated by energy gaps, as shown in the figure. As Npi become 
very large, the groups of states become continuous bands, with the photonic-lattice modes 
residing entirely within the shaded regions, and the free-space modes residing outside. At 
Tpi = 1, the model (with a very long L) approximates the free-space situation. 

Figure 4 shows that the model can also reproduce the significant fiattening of the 
photonic-lattice dispersion at the bandedges. Plotted in the figure is the dispersion for 
a twelve period {Npi — 12) photonic lattice, where the points indicate the actual eigenfre- 
quencies and the sohd curve is a fit of the data to illustrate the case of Npi — > oo. The 
fiattening of the dispersion at a bandedge results in a drastic increase in the density of states. 
We define the density of states as p{u) = dk^/dQ, where following solid state convention 
ko is the wavevector with vanishing interface refiectivity. The effects of the large density of 

states increase on intensity inside and outside of a photonic lattice is the focus of this paper. 
III. ACTIVE MEDIUM AND RADIATION FIELD 

To study the modification of emission characteristics by a photonic lattice, we consider 
the situation of an ensemble of two-level atoms located inside a photonic lattice. Each atom 
is labeled by n and j, so that |a„j) and \bnj) are the ground and excited states, respectively. 
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of an atom located at Zj inside the photonic lattice, with resonant energy fvjJn- Assuming the 
dipole approximation, Zj is a parameter locating the atom to a region that is small compared 
to a wavelength, but large compared to the size of an atom. We describe the radiation field 
emitted by these atoms in terms of the combined system eigenmodes derived in the previous 
section, i.e. 

E{z,t)^Y.Sk[ak{t) + al{t)\ Uk{z) (7) 

k 

where Sk — y/hflk/ (^^o), a^. and are the photon creation and annihilation operators, 
respectively, and A is the cross section area of the structure. Prom (7), using Maxwell's 
equations and a dipole interaction, the Hamiltonian for the matter and radiation-field system 

is20,21 

H = ^ huJn \b„j) {bnj\ + ^ Klka\ak - ^ Qkj (\hnj) (Onjl CLk + 4 lOnj) , (8) 

n,j k k,n,j 

where gkj = fiSkUk (%) and /i is the dipole matrix element. Introducing the operators for the 
microscopic polarization amplitude Pnjk = |6„j) (a„j| exp [—i {un — ^k) t], the excited and 
ground state populations, aanj = jo-nj) (flnjl and abnj = \bnj) {bnj\, respectively, and working 
in the Heisenberg picture,^* we derive the equations of motion 



dpnjk 


i 


dt 


^ h 




i 


dt 




dcTbnj 




dt 





-i{Wn—^k}t 



^fl'ifc'j (^CTbnjO'kO'k' ~ (A'^kCanj 



k 

5^ k 
njk 



) 



Pnjk^ 



Additionally, the photon number operator obeys. 



da\ak I 



dt 



t -i{uJj-nk)t _ i{u}j-Q.k)t 
Pnjk^ Pnjkc- 



(9) 
(10) 

(11) 



(12) 



Assuming that the polarization decays because of dephasing collisions and that the effective 
decay rate 7 is much larger than the rate of changes in the active medium and photon 
populations, we can adiabatically eliminate the polarization equation. Then, introducing 
the expectation values 
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N 



JV 



we obtain the working equations for our analysis: 



dNg, 
dt 



heoALcj 



Tfc [{Nbn - Nan) Nk + N^n] L (cU„ - Qfe) 



- 7r [-/Van - fa (c^n, 7")] - A {uj^i) K 



(it 



J2 Tifc [(A^6„ - Nan) Nk + A^6n] L {cOn - ^k) 



- Ir \Nlya - fb (oUn, T)] + A K) N, 

dNk 2/^2 



J2 ^k Tk [{Nbn - Nan) Nk + Ni^] L {uJn " ^k) " IcNk 



dt heoALc'f ^ 
where is the number of atoms, L (x) = [l + (2^/7)^] and 

/•Lc 

Tfe = / dz \uk {z)\^ 
Jo 



(13) 
(14) 

(15) 



(16) 

(17) 
(18) 



(19) 



is the mode confinement factor. In (16) - (18), the pump and decay contributions are included 
phenomenologically, 7c is the photon decay rate, A(a;„) = Aocxp [h{u!o — cUn) /ksT^ is the 
pump rate, hwQ is the material bandgap energy and 7^ is an effective rate for the actual 
populations Nan and A^n to relax to the equilibrium distributions 



fa{t^n, T) — Zn 
fb{uJn:T) = Z„exp 



— hWn 
1^ 



where 



1 + exp 



— hx)n 
1^ 



(20) 
(21) 

(22) 



Tp and T are the pump and reservoir temperatures. In our study, (16) to (18) are solved 
numerically. 



IV. DETECTOR 



To determine the spectra of the intracavity and output radiation, we use the simple 
spectrometer model shown in Figure 5. In this model, two- level atoms are placed in the 
region of interest. These atoms are prepared with only the ground state |a^) populated when 
the radiation field is absent (zero detector temperature) . The label n indicates that the level 
spacing between and |a^) is lo^. The atoms interact weakly with the radiation field 
to be measured, which excites some fraction of the atoms to an excited state \b'^) that has 
some finite lifetime 7^^. Assuming a sufficiently fast detector response so that the detector 
populations adiabatically follow the variations in the photon number, the population in state 
l^n) gives a measure of the radiation intensity (oc N^) in the region occupied by the detector 
atom. The steady state upper detector state population is 

where D — lii^Ndl ij^^-^^dld)-, fJ'd is the dipole matrix element between states and 
|a^), is the length of the detected region, and is the number of detector atoms. 
Measuring this population for atoms of different cu^ gives the spectrum within the region 
Zd ^ z < Zd + Ld- In this model, A^^ and the decay rate 75 should be sufficiently large to 
prevent saturation of the detector. On the other hand, too large a 7^ degrades spectral 
resolution. Alternately, one may use two level atoms injected into the region of interest, and 
removed after a short time.^^ 

V. PHOTONIC LATTICE EMISSION 

We consider a twelve-period photonic lattice with Lc = 120/im, L = 1.2cm and interface 
transmission Tpi = 0.01, 0.1 and 0.4. The eigenmodes are determined by solving (2) with 
the boundary conditions (3) - (5). The results are used in (16) - (18), which are solved 
numerically with a fourth-order Runge-Kutta finite difference method. The input parameters 
are 7 = lO^^s'^, -fc = 10^s-\ Aq = 10^°s-\ ujo = 1.6 x 10^^s-\ = e x 1.3nm, A^ = 601 
and Tp^T^ 400X. 
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A. Equilibrium 



To relate to earlier studies, we first compare photonic lattice and blackbody emis- 
sions under thermal equilibrium conditions. To do so, we perform the calculations for a 
rapid population relaxation rate of 7^ = lO^^s"^, which ensures (verified after the time 
integration) that the steady-state active-medium populations Nan and N^n are to a good 
approximation given by the equilibrium distributions fa{uJn,T) and fb{u!n,T), respectively. 
The solid curves in Fig. 6 (a) show the calculated intracavity emission spectra for three 
interface transmissions. In the figure, we define an intracavity detector signal, 

^ k Id + i^-^k) Jo 

where {cu) is calculated using the steady-state solution for A^^ in (23). The figure shows 
two bands of photonic-lattice states, where the frequency extent of the bands depends on 
the interface transmission. Between the two bands is a photonic bandgap where emission 
is strongly suppressed. By repeating the calculation with Tpi — 1, we obtain the corre- 
sponding blackbody spectrum (dashed curve). Comparison of the curves clearly indicates 
the significant intensity enhancement inside a photonic lattice, especially at the bandedges 
for Tpi = 0.01. 

To determine the output spectrum, we place the spectrometer in the free-space region. 
Figure 7 shows the output detector signal, 

for the same interface transmissions as in Fig. 6. In contrast to inside the photonic lattice, 
where there is significant optical intensity enhancement. Fig. 7 indicates that the intra- 
cavity emission peaks are appreciably depressed outside the photonic lattice. The strong 
intracavity enhancement by the photonic-lattice density of states appears to be cancelled by 
an outcoupling attenuation. This leaves the photonic lattice emission peaks to be essentially 
independent of interface transmission. More importantly, these peaks lie at or slightly below 
the blackbody emission curve. 

Simulations performed over a wide range of input parameters point towards the result 
that as long as the active-medium populations Nan and Nj,n are in thermal equilibrium, the 
photonic-lattice output is always below that of the blackbody. For instance, the spectra are 
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insensitive to the choice of 7. The blackbody spectrum is also insensitive to 7^ because of 
the weak frequency dependence of the blackbody photon density. However, for the photonic- 
lattice, too large a 7^ degrades the spectrometer resolution and leads to lower spectral peaks. 

For the opposite situation, too small a 7^^ introduces noise in the spectrum because of the 
inadequate resolution of the system normal modes (i.e. because L is insufficiently large). 

B. Nonequilibrium 

To study active photonic-lattice operation in greater generality, we allow the active 
medium populations to deviate from thermal equilibrium. The investigation is performed by 
repeating the earlier calculations, keeping all input parameters except 7,. the same. Figure 
8 illustrates the changes in the excited-state population distribution Nbn-, at steady state 
and for decreasing population relaxation. When 7^ is reduced to lO^^s"^ from 10^^s~\ a 
slight difference emerges between the excited-state populations of the identically pumped 
photonic-lattice and blackbody active media. The solid curve in Fig. 8 (a) shows a no- 
ticeable deformation of the photonic-lattice excited-state population distribution. There is 
also a significant difference between A^6„ and fb{uJn,T) for T = AOOK (dot-dashed curve), 
which are the actual distribution and the asymptotic (7^ 00) equilibrium distribution, 
respectively. Further reduction to 7^ = lO^^'s^^ significantly increases the deviation of the 
photonic-lattice excited-state population distribution from a Maxwell Boltzmann distribu- 
tion [see solid curve. Fig 8 (b)]. Holes are burned in the distribution because the population 
relaxation is insufficiently fast to replenish the excited-state population depleted by the spec- 
trally relatively narrow radiation field emitted by photonic lattice. There is also a change in 
the blackbody distribution [dashed curve. Fig. 8(b)], to one that approximates a Maxwell 
Boltzmann distribution at T ^ 500K (dotted curve). Since the active media in both struc- 
tures are identical, the difference between the photonic-lattice and blackbody populations 
(solid and dashed curves, respectively) is from photonic-lattice effects. 

The effects of the population changes in Fig. 8 on the emission spectra are depicted in 
Fig 9. Plotted on the y-axis is the relative emission intensity inside (outside) the photonic 
lattice, which we define as Sin{out) (^^) for the photonic lattice divided by S (uj) for the 
blackbody. In spite of the large increase in excited state population, we find that the 
intracavity and output relative intensities remain basically unchanged when 7^ is reduced 
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from lO^^s"^ to lO^^s"^. In particular, the output photonic-lattice intensity remains at 
or slightly below that of the blackbody (i.e., relative intensity < 1). However, the result 
changes considerably for 7^ — lO^'^s""^. Here, the intensity within the photonic- lattice band 
increases considerably relative to that of the blackbody both inside and outside the cavity 
(solid curves). More importantly, the solid curve in Fig. 9 (b) clearly shows greater output 
intensity for the photonic lattice than the blackbody throughout the emission band of the 
photonic lattice. This enhancement of output emission occurs for identically pumped active 
regions and is a result of a nonequilibrium population that shows significant hole burning. 
The presence of nonequilibrium effects may be the cause for experimental observations of 
metallic photonic-lattice emission exceeding that of the blackbody. Note that the difference 
in output emission spectra [solid and dashed curves in Fig. 9 (b)] comes from population 
distributions that are, on the average, quite similar. That is, a least-squares fit of the 
solid and dashed curves in Fig. 8 (b) will produce Maxwell-Boltzmann's distributions that 
differ in temperature by less than 20K. Therefore, measurement of average temperature 
will not identify the experimental conditions leading to the photonic-lattice output emission 
exceeding that of the blackbody. Rather, an energy-resolved measurement of the emitter 
upper state population is necessary. 

The end results reached in our analyses involving equilibrium and nonequilibrium sit- 
uations are robust, i.e., they are relatively insensitive to the choice of input parameters. 
While calculations performed with different interface transmission show significant differ- 
ences in spectral shapes and intracavity intensities, the output intensities remain relatively 
constant because of the mitigating influence of the coupling to free space. Calculations are 
also performed for different number of photonic-lattice periods. The results show negligible 
differences beyond N^i = 10, thus verifying that the use of a 12-period photonic lattice does 
not lead to loss of generality. Clearly noticeable are effects, such as differences in excess 
output intensity with varying interface transmission, that are due to optical nonlinearities 
in the nonequilibrium active medium. Such effects will not be present in treatments using 
linear classical sources. 
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VI. EXTENSION TO A 3-DIMENSIONAL PHOTONIC LATTICE 



This section treats a 3-dimensional photonic lattice. Following earlier quantum optical 
studies of photonics lattices, a spherically symmetric dielectric function is assumed to sim- 
plify the numerics. In spherical coordinates, the equation satisfied by the passive eigenmodes 
of the combined photonic- lattice and free-space system, Ukim {t^, d, <f>), is 

1 d^Uklr. 



d ( . ^dUklm 



+ 



sin^ 9 dcf)^ 



= -//oe (r) ^kim'^kim ■ (26) 



Choosing the dielectric function 

e (r) = eo 



(27) 



where rj and k are the same as in (1), a solution of (26) between the photonic- lattice interfaces 



or in the free-space region is 



Ukim (r, 9, (p) = [Akinji {kr) + Bkinm {kr)] Yi^ {9, (p) , (28) 

where ji (p) and rji (p) are spherical Bessel and Neumann functions, Yim {9, (p) is a spherical 
harmonic and the subscript n indicates that the coefficients Akin and Bkin are for r„ < r < 
Tn+i- In order for a solution to be finite at the origin and vanish at r = tn^i+i (the end of 
the region representing free space), we require 

Bkii = (29) 



AkiNj,i+iji (kr) + BkiNpi+iVi (kr) = . (30) 
At the interfaces, the boundary conditions (4) and (5) demand 

Akinji (kr) + BkinVi (kr) - Akm+iji {kr) - B^in+iVi (kr) = (31) 

Akin [-j'l (kr) + riJi (kr)] + B^n [-ti'i {kr) + rjrji {kr)] 

+Akin+ij'i {kr) + Bkin+iv'i {kr) = (32) 
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for 2 < n < Npi and r = ri, r2,...rN^i. 

The numerical solution is implemented similar to what is described in Appendix A. A 
6-period photonic lattice is considered, where the lattice constant is and the interface 
transmission is Tpi = 0.05. Coupled to the photonic lattice is a 'free'-space region extending 
from 24/im < r < 612/im. System dynamics is governed by the equations of motion (16) 
- (18) with the photon-state index k replaced by the three indices, k, I and m. Numerical 
analyses of the steady state solutions are performed assuming the input parameters, 7 = 
2 X 10^2^-1^ ^ 10^3s-^ 7c = 10^^s-\ Aq = lO^^s-^, uo = 1.6 x lO^'^s'^, // = e x 1.3nm, 
N — 601 and Tp = T = 200K. Similar to the 1-dimensional case, the rates are chosen 
to ensure reaching steady state with active-medium populations Nan and N^n described by 
Maxwcll-Boltmann distributions. 

To obtain the intracavity emission spectrum, detector atoms are placed inside the 
photonic-lattice structure. The probability of finding a photon with frequency cu is pro- 
portional to 

Sin{uj) ^\^Q.kimNkim—^ ^9 1 dr r"^ \ukim {r, , , (33) 

tiZ Id + i^- ^klm) Jo Jo Jo 

so that the emission energy is proportional to hjjjSin{oj). In Fig. 10, the solid curve is a 
plot of ujSin as a function of frequency. It shows four narrow emission bands separated by 
photonic bandgaps. Repeating the calculation using interface transmission Tpi — 1, gives 

the free-space emission spectrum (dashed curve). Examination of the populations after 
steady state is reached verifies that both solid and dashed curves arc for identical Maxwell- 
Boltmann distributions at T = 200i^. The curves clearly indicate intensity enhancement 
inside the photonic lattice. 

To compare the output emission in a given direction, the detector atoms are placed to 
give a signal. 



'S'o«t(u;) = V'^^WmA^feim^ , ^'^ — 72 I d<P I dO i dr r"^ \ukim{r,e,(j))\^ , 

klm 

^''klmj Jtpdi "^^di ^"i^ 

where r^i, rd2 are within the free-space region, 0,^, 9^2 and 4>di, 4>d2 define the direction and 
collection sohd angle. The sohd curve and dashed curves in Fig. 11 shows the output 
photon-lattice and free-space emission spectra, respectively. These results are obtained for 
rdi — lOO/im, rd2 — 124//m, = 9di — 0, 0^2 = 27r and 9d2 — 7r/18, which define emission 
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d2 



rd2 



within a cone of ±10° in the z-direction. Comparison of sohd and dashed curves indicates 
that, similar to the 1-dimensional case, peak intensities measured outside the photonic-lattice 
structure do not exceed those of the blackbody. 

A test of our treatment is to see how well it reproduces Planck's distribution for the 
frequency spectrum of free-space emission energy from a thermal source. The dotted curve 
in Fig. 11 is proportional to uj"^ [exp (hw/kBT) — with temperature T = 200K, and 
it depicts the shape of the blackbody frequency spectrum according to Planck's formula. 
The agreement is good, considering that there are several factors causing discrepancies. 
Two important ones are truncating the optical modes at Z = 13 and limiting 'free'-space to 
24//m < r < 612/xm, in order to maintain reasonable computation times. Even so, over 10^ 
optical modes are used. Other factors contributing to the differences include the presence 
of optical loss [7c 7^ in (18)], which is neglected in the derivation of Planck's distribution. 
While increasing the number of optical modes improves agreement with Planck's formula, 
it does not impact the blackbody versus photonic lattice emission comparison. This is 
because the result of the photonic-lattice output intensity spectrum being bounded inside 
the blackbody one applies separately for each I. 

VII. VI. CONCLUSION 

In summary, the emission from an active photonic lattice is investigated using a model 
consisting of an inhomogeneously broadened ensemble of two-level atoms interacting with a 
multimodc radiation field. A fully quantized (i.e., quantized atoms and quantized electro- 
magnetic field) description is chosen to provide a consistent description of stimulated and 
spontaneous emission. Furthermore, to describe the modal properties of the radiation field 
of a finite photonic lattice coupled to free space, the analysis considers the photonic lattice 
and free space as one combined system. This circumvents a long-standing inconsistency in 
quantum optics involving the decoupling of the treatments of the cavity normal modes and 
out coupling losses. 

Our approach gives the emission spectra for arbitrary photonic-lattice configurations 
and reproduces Planck's blackbody radiation formula for thermal emission in free space. 
Comparison of photonic-lattice and blackbody emission shows appreciable modification of 
the blackbody spectrum by the photonic lattice, where the redistribution of the photon 
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density of states results in suppression of radiation at certain wavelengths and enhancement 
at others. The enhancement can give rise to high intracavity intensity peaks, especially at 
the photonic-lattice bandedges. These intensity peaks are mitigated outside the photonic 
lattice by the spectrally dependent outcoupling. For population relaxation sufficiently fast to 
ensure the same equilibrium population distribution in both structures, the photonic-lattice 
output intensity does not exceed that of the blackbody at the same frequency. However, 
for slow population relaxation, there is a greater tendency for a nonequilibrium photonic- 
lattice population. Then, in the presence of population hole burning, the intensity in certain 
regions of the photonic-lattice spectrum can exceed that of the blackbody, even when both 
structures arc identically pumped. 
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IX. APPENDIX A: NUMERICAL EVALUATION OF SYSTEM EIGENMODES 

In this appendix we describe a numerical procedure for evaluating the eigenmodes of 
the combined photonic-lattice and free-space system. This procedure applies for a photonic 
lattice of arbitrary size and interface transmission. In region n, which may be any section 
of the photonic lattice or the section representing free space, the solutions of (2) have the 
form 




(Al) 



where k — Vlk/c- Because of boundary condition (3), 



Ak,i ^ 



(A2) 



(A3) 



^fc,jv^,+i sin [kzN^i+i] +B, 



(A4) 
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Prom boundary conditions (4) and (5), 



A, 



.fe,n sin {kzn) + Bk,n COS {kzn) = ^fe,n+i sin {kZn) + B, 




^k,n [- COS {kZn) + sin {kZn)] + -Bfc,„ [sin (A;^;^) + 77 cos (A;^;^)] = -A^^n+i cos + -Bfc,n+i sin {kz^) 



for 2 < n < A^pi. 

There are many approaches to numerically solve the above equations. We describe below 
the one we followed in this paper. Basically, we look for the values of k satisfying (A4), 
where the mode amplitudes AA;,Arp,+i and -Bfc,iVpi+i are obtained by solving a 2Npi x 2Npi 
matrix equation 



The matrix elements of S are best defined by separating the even and odd number rows. 

For n = odd, 



(A6) 



SU 



(A7) 



Si^i = - sin (kzi) 
Si,i+i = - cos (kzi) 



(A8) 



(A9) 



and for n = even, 



sin (kzi-i) 



(AlO) 



For n > 3 and n — odd. 



= cos (kzi) 
= sin (kzi) 



(All) 



(A12) 



and for n > 4 and n — even. 



Si^i_i = cos {kzi_i) 

= sin (kzi-i) + 7] cos (kzi-i) 
Si,i-3 = - cos (kzi^i) + 77 sin (kzi^i) 



(A13) 
(A14) 



(A13) 
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All other matrix elements are zero. The elements of the column matrix D vanishes except 
for 

Di = - sin {kzi) (A15) 
D2 — cos{kz2) — r]sva.{kz2) (A16) 

Equation (A7) is solved using the Gauss- Jordan method. In the solutions and for j — odd, 
Uj gives the coefficient Afcj+i, while Uj+i gives the coefficient Bkj+i. At this stage, we have 
set 1 = 1 so that the eigenfunctions are unnormalized. We perform the normalization 
according to (6). 
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Figure Captions 



Fig. 1. Model of a photonic lattice connected to a large cavity approximating the universe. 

Fig. 2. Eigenfunctions of a 6-period photonic- lattice coupled to free space for interface 
transmission Tpi — 0.10. The figure shows non-resonant (a) and resonant (b and c) photonic- 
lattice modes. 

Fig. 3. Eigenfrequency versus interface transmission for a 6-period photonic lattice. 
The points indicate the actual cigcnmodes of the finite structure, while the shaded regions 
illustrate the extent of the bands of photonic-lattice states when the number of periods 
become very large. 

Fig. 4. Dispersion for a 12-period photonic lattice. The points are the eigenmodes, the 
solid curve is a fit through these points and the dashed curve shows the free-space dispersion. 

Fig. 5. Spectrometer model where the upper level decay prevents detector saturation 
and approximates the drift of carriers to the electrodes in a reverse-baised photodiode. 

Fig. 6. Photonic- lattice (solid curves) and blackbody (dashed curves) intracavity emission 
spectra. 

Fig. 7. Photonic-lattice (sohd curves) and blackbody (dashed curves) output emission 
spectra. 

Fig. 8. Upper state population for photonic lattice (solid curve) and blackbody (dashed 
curve) versus transition frequency for 7^ = 10^^ (a) and 10^*^5'^ (b). The dot-dashed curves 
show the equilibrium distribution at 400i^. The dotted curve is the equilibrium distribution 
at 500K. 

Fig. 9. Relative intensity spectra inside (a) and outside (b) the photonic lattice in Fig. 
8. The curves are for 7,. — 10^^ (dashed) and 10^°s~^ (sohd). Above the long-dashed line, 
the photonic-lattice intensity is higher than the blackbody's. 

Fig. 10. 3-dimensional photonic- lattice and blackbody intracavity emission spectra (solid 
and dashed curves, respectively). 

Fig. 11. Output emission spectra for 3-d photonic-lattice (solid curve), blackbody (dashed 
curve) and Planck's distribution (dotted curve). 
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